Probing local relaxation of cold atoms in optical superlattices 
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In the study of relaxation processes in coherent non-equilibrium dynamics of quenched quantum systems, 
ultracold atoms in optical superlattices with periodicity 2 provide a very fruitful test ground. In this work, we 
consider the dynamics of a particular, experimentally accessible initial state prepared in a superlattice structure 
evolving under a Bose-Hubbard Hamiltonian in the entire range of interaction strengths, further investigating 
the issues raised in Ref. [Phys. Rev. Lett. 101, 063001 (2008)]. We investigate the relaxation dynamics ana- 
lytically in the non interacting and hard core bosonic limits, deriving explicit expressions for the dynamics of 
certain correlation functions, and numerically for finite interaction strengths using the time-dependent density- 
matrix renormalization (t-DMRG) approach. We can identify signatures of local relaxation that can be accessed 
experimentally with present technology. While the global system preserves the information about the initial 
condition, locally the system relaxes to the state having maximum entropy respecting the constraints of the ini- 
tial condition. For finite interaction strengths and finite times, the relaxation dynamics contains signatures of 
the relaxation dynamics of both the non-interacting and hard core bosonic limits. 



I. INTRODUCTION 

Both in classical and quantum physics, a complete frame- 
work for the description of equilibrium properties of arbitrary 
physical many-body systems exists, although the explicit cal- 
culation of many-body equilibrium properties is an often un- 
solved challenge. The situation is much less satisfactory when 
it comes to the study of the non-equilibrium properties of 
many-body systems, where such a general framework is miss- 
ing and may not even exist. Research in this field has there- 
fore focussed on relatively specific issues and types of non- 
equilibrium. One of the issues taking center stage is whether 
quantum many-body systems in non-equilibrium evolving co- 
herently under a local Hamiltonian equilibrate or not. If so, 
one may ask whether the equilibrium states can in some sense 
be described by a thermal state. This old and fundamental 
question of the equilibration of quantum many-body systems 
has enjoyed quite a renaissance recently [1-36]. 

A specific setting of coherent non-equilibrium quantum dy- 
namics is provided by quantum quenches, where one starts 
from an eigenstate of some initial Hamiltonian and pushes the 
system out of equilibrium by a sudden change (or "quench") 
of system parameters, leading to a new Hamiltonian. One then 
considers the evolution of the system under the new Hamilto- 
nian. A further restriction is provided by the assumption that 
the quantum system under consideration is closed, i.e., has no 
coupling to a bath of degrees of freedom that might assist the 
relaxation process. Time evolution (and hence the potential 
relaxation to some equilibrium) will obviously be constrained 
by the constants of motion, i.e., Hermitian operators commut- 
ing with the new Hamiltonian whose expectation values are 
fixed by the initial state; in that sense, any relaxed state will to 
a certain degree show some memory of the initial state. 

It has been conjectured that — in some sense — the quantum 
system should relax to the maximum entropy state consistent 
with the expectation values of the constants of motion fixed 
by the initial state (see, e.g., Refs. lfT8lfT9ll27l ). also referred 
to as a generalized Gibbs ensemble fyT\. This may be attrac- 



tive because it reminds of Jaynes' derivation of equilibrium 
statistical mechanics. 

This observation is in conflict with the fact that obviously, 
if the system can be meaningfully treated as a closed quantum 
system, one cannot expect the whole system to relax: Initially 
pure states will remain so in time under a unitary time evolu- 
tion [2-6]. After all, the entire information about the initial 
condition is still stored in the system, albeit in a dilute fash- 
ion. Yet, this observation is by no means in contradiction with 
the possibility that in any local observation, the system may 
appear perfectly relaxed, even without invoking a time aver- 
age 12] El El |6l . The key point is that locally, one may well 
expect the relaxation to be true [38] : For any subset of sites in 
a sufficiently large lattice system, the reduced state may well 
converge to the reduced state of the maximum entropy state, 
given the conserved quantities of motion, and stay relaxed for 
an arbitrary long time. Indeed, for such a subset it is, under 
suitable assumptions about its interactions with the rest of the 
world, easy to make contact to Jaynes' formulation of statisti- 
cal mechanics. 

Ref. [3] considers a variant of the question in which this 
local relaxation of subsystems, referred to as local relaxation 
conjecture, can in fact be rigorously proven to hold exactly: 
This is the one where one evolves a state deep in a Mott phase 
according to a Bose-Hubbard Hamiltonian coiTesponding to 
the deep superfluid phase — treated as a non-interacting sys- 
tem. In this setting, the local relaxation is in fact true: The re- 
duced state of a block / of consecutive sites indeed converges 
to a maximum entropy state 

(1) 

(in trace norm) for large systems and large times, having max- 
imum entropy consistent with the constants of motion [3|. 
Note that there is no time average and the initial state is not a 
Gaussian state. 

Also, for free bosonic and fermionic systems, and for Gaus- 
sian initial states, it has been shown rigorously in Refs. ||5l|6l 
under which conditions the local relaxation conjecture is true. 



It turns out that in these cases, the presence and speed of local 
relaxation depends crucially on the single-particle spectrum 
and the dimensionality of the systems. 

The in some sense inverse case of a quench from the super- 
fluid phase to the Mott insulating phase, but generally at finite 
interaction strengths, has been studied in Ref. 121. In this non- 
integrable system, numerically two distinct non-equilibrium 
regimes have been found where equilibrated states resembled 
thermal states or states with memory. 

The physical intuition why local relaxation happens is the 
following: If one switches to a new Hamiltonian, the system 
is no longer in equilibrium. Hence, one has local excitations 
at each point [2, 3, 12-16]. They cannot travel arbitrarily 
fast, however, as there is a finite speed of information trans- 
fer in lattice systems 1391 l40l . At each site, in time more and 
more "waves" of farther and farther away sites can possibly 
have a significant influence on this site. This is related to a 
finite speed of sound or of information transfer in a lattice 
system |39|. Then, local relaxation may be a consequence of 
the incommensurate influence of these excitations generating 
mixing and thermalization. The "thermalization" time scale is 
hence governed by the speed of information transfer |46|. 

This also links to kinematical approaches to the problem 
[37, 47-51], arguing that most states anyway look locally very 
much relaxed in that they have large entropy. A random pure 
state (as taken from the Haar measure) will locally have a large 
entropy. Specifically, in interacting systems one should expect 
such a local relaxation to be true, too, an aspect that will be 
studied in this work. 

So far the discussion has not taken into account the recur- 
rence happening in any closed quantum or classical system 
131. For finite, but large system sizes, recurrence times will 
however become so long that recurrence effects become neg- 
ligible, in that the quantum many-body system can be effec- 
tively locally equilibrated on a much shorter time scale than 
the recurrence time. 

We will see that while settings exhibiting local relaxation 
may be generated in various fashions, it is a greater challenge 
to actually probe signatures of local relaxation. This apparent 
dilemma — that to demonstrate local relaxation appears to ne- 
cessitate local addressing — will be resolved in this work, by 
making use of a period-2 setting, further developing the idea 
of Ref. |4 |. The taken path opens up a way to experimentally 
explore local relaxation effects using atoms in optical super- 
lattices. We systematically investigate various aspects of local 
relaxation in such a setting, both using analytical as well as 
numerical methods, based on a time-dependent density-matrix 
renormalization-group (t-DMRG) approach. 



II. EXPERIMENTAL SETUP: ULTRACOLD ATOMS IN 
OPTICAL SUPERLATTICES 

The current surge of interest in the relaxation of quantum 
systems after a quench is mainly motivated by the advent of 
ultracold atoms in optical lattices [52-56]. These systems are 
highly attractive as they allow for sudden controlled manip- 
ulations of system parameters, hence quenches, are strongly 
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FIG. 1. Idealized sketch of the effect of local relaxation in a setting 
having periodicity 2. Starting from an idealized initial condition of 
bosonic atoms being present or absent, |i/'(0)) ~ |1, 0, . . . , 1, 0), 
in even and odd sites of a one-dimensional Bose-Hubbard system- 
achieved by imposing a superlattice to the one-dimensional optical 
lattice-the system locally relaxes to an apparent maximum entropy 
state. 

interacting, hence non-trivial, are on experimental time scales 
essentially closed quantum systems, and therefore show co- 
herent quantum dynamics. Systems are also sufficiently large 
to show non-tiivial many-body behavior. 

In the present context, however, the major drawback of ul- 
tracold atoms is that despite the unprecedented possibilities of 
manipulation the study of local relaxation provides an exper- 
imental challenge. This is due to the fact that local (i.e., site- 
resolved) measurements on ultracold atoms in optical lattices 
are still not satisfactory, albeit rapid progress is being made. 
In this work, we propose to study instances of local relaxation 
even in a setting where strictly local quantities can not easily 
be studied by using optical superlattices [57-60]. 

Following the setup very recently realized by Bloch and 
coworkers 157 . 58] , we consider bosonic Rubidium-87 atoms 
in a period-2 optical superlattice geometry: Two standing- 
wave laser fields at wavelengths 765nm and 1530nm are su- 
perimposed with fixed relative location to provide a superlat- 
tice geometry of lattice constants a = 382. 5nm and 2a = 
765nm. It is experimentally possible, among other things, to 

(i) change the relative strength of the two optical lattices and 

(ii) shift their relative position, by altering phase and detuning 
of the optical superlattice. Assuming the strength of the lat- 
tice with lattice constant a to be fixed, (i) allows to couple and 
decouple this lattice into an aiTay of double well potentials, 
combining one odd (o) and one even (e) site of the original 
lattice. In such double well potentials, (ii) allows to introduce 
a bias between the chemical potentials of the odd and even 
sites, A = — /ig. 

Isolating double wells and biasing odd vs. even sites, in 
turn, allows for the preparation of patterns of atoms and for 
extracting local quantities to the degree that odd and even sites 
can be distinguished from each other: 

• Preparation of periodic patterns is achieved by loading 
the superlattice while introducing a bias between odd 
and even sites such that due to a shifting of particles 
all particles are on either odd or even sites after load- 
ing. Using further experimental techniques [57], multi- 
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pie occupancies on the occupied sites can be eliminated, 
leaving a sequence of empty and single-occupied sites. 

• Period-2 local measurements can be obtained by map- 
ping odd and even sites to different Brillouin zones: As- 
suming completely decoupled double wells, each part 
of the double well has multiple bands separated by well- 
defined energies. Biasing, say, the odd sites relative 
to the even sites by an energy in excess of the separa- 
tion energy of the band-separation energy, odd-site par- 
ticles are reloaded into the higher band of the even sites, 
whereas the even-site particles stay in the lower band. 
A standard time-of-flight mapping then shows the even- 
site particles in the first Brillouin zone, the odd-site par- 
ticles in the higher Brillouin zones. 

• More sophisticated measurements with period-2 can 
also be performed in principle. When letting the sys- 
tem evolve for some defined hold time, and then freez- 
ing the time evolution by ramping up the barrier, one 
can also measure nearest-neighbor correlations in the 
lattice |B71 - |601 . There are several ways to proceed here: 
On the one hand, one can isolate double wells, let them 
dephase. Then upon free expansion, the average corre- 
lations can be measured as in a double slit experiment. 
One the other hand, correlations can be mapped to den- 
sities. To the extent that the barrier between the wells 
is sufficiently high such that the time evolution can be 
described by a collection of double well systems (and 
higher Wannier bands do not have to be taken into ac- 
count), we can investigate the time evolution in each 
individual double well. Up to the presence of a confin- 
ing potential, the total time evolution then reflects the 
time evolution in each of the double wells. Let us label 
the two modes in any of the wells by 1 and 2, one may 
apply an appropriate free Hamiltonian to map correla- 
tions onto on-site densities. Specifically, for L/ = we 
find that 

2im(ala2> = (4Wa2(t)) - {a\{t)di{t))), (2) 
when letting di, 0,2 evolve under the free Hamiltonian 



h = — (^d\d2 + aldi ) , 



(3) 



until t — 7r/4. One can hence measure the imagi- 
nary part of the appropriate correlators, and hence the 
period-2 correlators in the full lattice. In settings where 
in the Heisenberg picture the phase map 0,2 1— >■ 102 for 
each right well is feasible, one can also measure the 
real part. If one has small interactions (because of the 
Wannier band problem — a description in terms of Wan- 
nier functions still has to be valid, accompanied by a 
non-vanishing U), then this leads to a small dephas- 
ing in this mapping. Direct numerical simulation can 
take this properly into account, allowing for a correct 
interpretation of the experimental observation. Such a 
technique allows to measure correlators from a time- 
resolved observation of on-site densities. Similar time- 
resolved measurements have recently been performed 
in optical superlattices Ii58ll . 



As shown in Fig.[T| we propose to start from a two-periodic 
initial state prepared by the superlattice setup as described 
above, where all odd sites are occupied by exactly 1 boson 
and all even sites are empty. The initial state vector of the 
entire lattice is hence 



1^(0)) = |1,0,1,0,...1,0). 



(4) 



Tools how to experimentally achieve that sites are to a good 
approximation occupied by a single atom, and not two or more 
atoms, are described in Refs. [57-59]. If the 2a-lattice is sud- 
denly switched off, i.e., the system quenched, the state vector 

conventional Bose-Hubbard Hamiltonian 



will evolve in time \i^{t)) = c ^^\ip{0)) according to the 



H 



U 



^nj(ni-l), (5) 



where U and J are the standard interaction and hopping pa- 
rameters of the Bose-Hubbard model that can be calculated 
microscopically from the lattice parameters. The system size 
is given by L, whereas the boson number is = L/2 from 
this setup. L will always be even in the following, in line 
with the proposed setup; various boundary conditions will be 
imposed. We work in units where J = 1 and h= 1. 

We will not consider the effect of the occupation of 
higher order bands, but will stay within the limit of appli- 
cability of the Bose-Hubbard model. We will for the pur- 
poses of the present paper also neither consider an additional 
harmonic confining potential — except when explicitly stated 
otherwise — nor additional dephasing effects due to statistical 
fluctuations of local fields. Instead, we will systematically 
flesh out what behavior is expected in this slightly idealized 
setting, to see how local relaxation manifests itself here, to 
comment on the impact of imperfections later. 

As ^(0) = |'0(O))('0(O)| is not an eigenstate of H, we are 
in a non-equilibrium situation. The initial state is uncorrected 
and shows inhomogeneous density of periodicity 2. Coherent 
quantum dynamics is now expected to homogenize densities 
locally and to build up non-local correlations. 

Local relaxation is now monitored by the global mea- 
surement of the total occupation of even, (7Ve(<)), and odd, 
{No{t)), sites. In a translationally invariant setting, this gives 
access to local observables as 



L/2 



mt)) = Y,{h2dt)) = ^{n2^{t)), 



(6) 



and similarly for odd sites. 

If one has experimental access to the variances a^^^ e o ^ 

i-^tot e o) ~ {^tot,e,o)'^, dcnsity-density correlations (see also 
Ref ll6ll ) between all even and all odd sites may be obtained 
through 

2 

L/2 

= {{f^2^{t)n2J-l{t)) - {h2,{t)){h2j-im . (7) 

i,j = l 
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In fact, the value of Ntot is fixed to L/2 for the proposed pat- 
tern, so in repeated experiments of same length, this quantity 
will have variance zero, a^g^ — 0. 

Moreover, the quasi-momentum distribution defined as 



^(^'O = Z £ "^'^^'^ ' 9 e [0,27r], (8) 

is also accessible as a global quantity, from time-of-flight 
measurements. 

We now proceed as follows. In order to obtain analytical 
results we set [/ = and U = oo respectively, which are both 
essentially non-interacting limiting cases, can be solved ex- 
actly and show very similar, but not generally identical results 
(compare also Ref. il3l ). These results will be extended to the 
interacting finite-C/ case which will be studied using t-DMRG 
method. For the time-scales considered it is an effectively 
quasi-exact method. 



III. EXACT SOLUTIONS FOR LIMITING CASES 

Both limiting cases U = and U ~ oo are or can be 
mapped to free models. This case has been considered for 
Gaussian initial conditions in Refs. |5, 6|, as well as for non- 
Gaussian product initial conditions in Ref. ID. In these cases, 
the reduced state of a subsystem / 



(9) 



will converge for large systems and long times to ^max. be- 
ing the maximum entropy density operator as constrained 
by the initial conditions (or, more precisely, if recurrences 
are considered, this will be true in 1-norm to an arbitrarily 
small approximation error for an arbitrarily long time). Also, 
fermionic Gaussian initial states have been considered in Refs. 
|l3]|6l. Local relaxation will therefore happen in these cases. 

The physical mechanisms behind the relaxation processes 
in the cases ofU = Q and [/ = oo is sightly different, how- 
ever: In the case U = 0, it is due to dephasing in the sense 
that freely propagating bosons lead to reduced density oper- 
ator contributions of quickly oscillating phases that average 
out. In Ref |3|, this intuition leads even for non-Gaussian ini- 
tial states to maximum entropy states, by invoking a quantum 
version of a central limit theorem, and exploiting the finite 
speed of information transfer 1 39 1 . 

In the case U = oo, real scattering processes happen, albeit 
of a very specific form that allows for a formal mapping to 
a non-interacting fermionic problem (and to the XX model, 
compare also Ref. 1 13 1): The scattering is simply relegated to 
the internal sign structure of the wave function. We will also 
see that the two relaxation processes lead to different time- 
evolutions of most physical observables. 

For all cases of nonzero finite U, the situation should be 
somehow intermediate. For the specific setup chosen here, 
interacting particles will learn of each others existence only 
after some initial time they need to come into contact. We 
would therefore expect that for very short time scales, observ- 
ables should evolve as in the U = limit, with a crossover 



in behavior for longer time scales when they start interacting. 
The question is whether for which interaction strengths the 
limiting pictures remain essentially valid and whether there 
is a genuinely different intermediate interaction regime with 
different relaxation behavior. Compared to Refs. |3 5 6|, in 
this work we focus to a lesser extent on rigorous mathemati- 
cal convergence statements — like invoking quantum versions 
of central limit theorems — but instead put more emphasis on 
the phenomenology of the physical relaxation process as such 
in the 2-periodic setting. 



A. Non-interacting bosons: U = 

In the translationally invariant case of non-interacting 
bosons, the Hamiltonian can be diagonalized by Fourier trans- 
forming to new bosonic operators ai , . . . , to yield 



ak, 



k=l 



where 



A. 



-2cos(27rA:/i) 



(10) 



(11) 



are the eigenvalues of the circulant Hamilton matrix H with 
entries Hjj = — 5dist(ij),i- While the real experiment will 
not have periodic boundary conditions, our calculations show 
that for realistic system sizes and time scales the difference 
between open and periodic boundary conditions is negligible. 

In the Heisenberg picture, the initial state vector |^/'(0)) = 
|1, 0, . . . , 1, 0) remains time-independent, whereas the opera- 
tors evolve in time as 



kit) 



1 ^ 
-T 



(12) 



g-itAfcg27rifc(j-j)/L 



Straightforward algebra then yields the exact time-evolution 
of two-point correlations, see Appendix, 



(13) 



In the thermodynamic limit L oo this turns into 

(-l)'iJ-* 



Ik, 



-Jj-,{Ut), 



(14) 



where J„ is a Bessel function of the first kind. 

From the above expression, one can derive the quasi- 
momentum distribution for finite L, see Appendix. For 
Lql{2-K) e {1, ... ,L} it is constant, S{q,t) = 1/2, while 
one finds 



S{q,t) 



1 



L 

p=i 



,4Jitcos(27rp/L) 

sin(27rp/L — q) 



(15) 
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for all other q€ [0, 27r]. 

With slightly more effort, density-density correlations as in 
Eq. (|7]i emerge as (see Appendix) 

{n2^{t)n2j-i{t)) - {n2^{t)){h2j-i{t)) (16) 

L/2 

— /2i,2j-l/2j-l,2i — |V2i,2fc-l(^)V2j-l,2fc-l(i)P, 
fc=l 

which, as they are local quantities, relax for large systems 
and long times. For large L one finds for the global density- 
density correlator (see Appendix) 

{Ne{t)No{t)) - {Ne{t)){No{t)) 

L , „ (17) 

^ "Y^ (3 + Jo(8Ji) - 4[Jo(4Jt)]2) . 

We also show in detail the exact local relaxation to a 
maximum entropy state for the case of the initial condition 
1^(0)) = |1, 0, . . . , 1, 0), largely following Ref. H: For ev- 
ery block of consecutive sites / = {l,...,s}, every small 
error e > and any desired, arbitrarily long recurrence time 
tree > there exists a system size L and a local relaxation 
time ti-ei > such that 



\\Ql{t) - Qn 



< e 



(18) 



for all times t e [ti-ei, irec], see Appendix. Here, || • ||i denotes 
the trace norm. Hence, after the quench for U — for 2- 
periodic initial conditions, the local state becomes in time a 
maximum entropy state, to an arbitrarily good approximation. 



B. Hardcore bosons: U — > cxa 

In the limit (7 — > cxo, the interaction manifests itself exclu- 
sively in that bosonic occupation numbers are upper bounded 
by 1. This leads to a well-known mapping in case of quan- 
tum ground states: The hard core limit is equivalent to the 
XX spin model and a model of spinless free fermions. But 
even in time evolution, in the limit of large U, the population 
of sites with particle number larger than unity is dynamically 
suppressed to an arbitrary extent: The expectation value of the 
new Hamiltonian is obviously preserved under time evolution, 
tr[i?^?(0)] = ix[HQ{t)], H = -Ho + UHi, where 

L ^ L 

Ho = Y.^b\+ih + blk+i), ^1 = 2 E - 1) ■ (19) 

i=l i=l 

For initial states that are supported on span(|0), |1)) one 
has t.T[Hg{0)] -iv[Hog{0)]. Furthermore t.T[Hg{t)] > 
tr[{~Ho + UH2)Q{t)], where 



L oo 



H2 = Y.Y.\^){ki 



(20) 



i=l k=2 



where the latter projector defined on each site i. This, in turn, 
means that 



tr[go(e(0 - g(0))] 

u 



> \i[H2git)] > 0. 



(21) 



Hence, one can ensure to arbitrary accuracy for large U that 
g{t) is only supported on span(|0), |1)). One therefore ar- 
rives in a perfectly meaningful way at the familiar hard core 
limit of the Bose-Hubbard model, which is equivalent to a free 
fermionic spinless system. 

This mapping to non-interacting spinless fermions is done 
through the familiar Jordan- Wigner transformation 



(22) 



Under this transformation the initial state vectors turn into 

IV) =cl(0)cj(0)c^(0)...|0) (23) 

and the Hamiltonian reads upon a Fourier transformation to 
new fermionic operators di, . . . ,(1^ 



Hu=oo — E ^kd\dk 



(24) 



k=l 



with Afc as in Eq. (Ill and the time-evolution of the operators 
is given by 



(25) 



with Vi,j{t) as in Eq. (12 1 



This is formally identical to the 
U = Q case. However, there are two differences: Peri- 
odic boundary conditions in the original bosonic model map 



to (anti-)periodic boundary conditions ffj^j^^ = (^1) 
Hence, boundary conditions stay periodic if the particle num- 
ber N is odd, to which we restrict ourselves in the following. 

The second, more important difference is that differences in 
results will show up due to the Jordan- Wigner transformation 
of operators and the difference between bosonic and fermionic 
(anti)commutators. It is easily shown that local densities and 
correlations between nearest neighbors translate directly (see 



XIII D I, such that the respective results for [/ = oo are iden- 
tical to those for [/ = 0. For longer-ranged two-point correla- 
tions, including structure functions, as well as density-density 
correlations results differ. Density-density correlations are 
given by (see |XIIID] i 



{h2i{t)n2]~i{t)) - {n2i{t)){n2j-i{t)) 



/2i,2i-l/27'- 



(26) 



2j-l,2i, 



where fi j is as in Eq. ( 13 i. For large L the global density- 



density correlator is then given by (see Appendix) 
^-^(1- Jo(8Ji)). 



(27) 



C. Interacting softcore bosons: Finite U 

In this case, we are no longer facing an integrable model. In 
order to study the relaxation dynamics, we will therefore turn 



to the time-dependent variant (64^-'66] of the DMRG (density- 
matrix renormaUzation group) method l62ll63l . This method 
allows to follow the coherent time-evolution of strongly in- 
teracting quantum systems very precisely. Its reach in time 
is, however, limited by the growth of entanglement in quan- 
tum systems: Linear entanglement growth, for example, leads 
to an exponential growth in numerical resources needed. As 
we will see, interestingly, the system under study is charac- 
terized by very strong linear entanglement growth. In the 
free instances this linear increase is indeed provably correct 



see Sec. IX Hence, the reachable times are quite short 



(Jt ~ 6) with up to 5000 states kept in the simulations. For 
most quantities, this turns out to be sufficiently long to make 
contact to the non-interacting results and to read off long-time 
behavior. In particular, the results allow for a quantitative 
comparison to experiments. 

A subtlety arises from the fact that DMRG prefers open 
boundary conditions, and is hence closer to experiment. How- 
ever, we would like also to compare to analytical results, 
where periodic boundary conditions are preferable. As it will 
turn out, for system sizes and times considered, the difference 
is negligible. 



IV. TIME EVOLUTION OF DENSITIES 

For the time evolution of local densities both exactly solv- 
able cases U = and U — oo give the same result 



1 (-1)' 

2 2L 

1 (-1)' 



k=l 



,4itJcos(2Trfc/L) 



(28) 



Odd- and even-site densities relax symmetrically about the 
n = 1/2 axis to n = 1/2, with an asymptotic decay as 

t-l/2+o(t-l/2). 

All <-DMRG results for finite U are perfectly compatible 
with relaxation of local densities to n = 1/2. On very short 
time scales {t < 1) particles have typically not collided yet 
and are not yet sensitive to the values of U, hence, finite U re- 
sults are identical to the limiting cases [/ = 0, oo. The relax- 
ation behavior for intermediate times, however, deviates quite 
strongly. Interaction effects become visible right after parti- 
cles make contact. This can be seen in Figs. [2]|6]( all calculated 
for L — 32), where we compare the time evolution of local 
densities for various finite values of U to the special cases 
U = 0, oo. For small U (exemplified hy U = 1.5) and large 
U (exemplified hy U — 8), the comparison indicates that the 
relaxation of local densities is governed by the behavior of the 
limiting non-interacting cases. For intermediate values of U 
(exemplified hy U = 3), scattering seems to be most effective 
and lead to much faster damping and relaxation, much beyond 
the above square root time dependence. This is plausible, as 
close to a quantum critical model there is a limiting point in 
the spectrum of the new Hamiltonian, leading to specifically 
effective relaxation. Deviations from the limiting behavior are 
sufficiently strong that they should be visible experimentally. 



U=1.5 (even) 
U=0 (even) 
U=1 .5 (odd) 
U=0 (odd) 




FIG. 2. Local density {ni{t)) vs. time, showing local relaxation. 
Shown is the time evolution of an even and an odd site for [/ = 
and U = 1.5. 



U=3 (even) 
U=0 (even) 
U=3 (odd) 
U=0 (odd) 




FIG. 3. Local density {hi{t)} vs. time, showing local relaxation. 
Shown is the time evolution of an even and an odd site for [/ = 
and U — 3. Note the strong deviation from the non-interacting limit 
and the strong suppression of density oscillations. 



U=8 (even) 
U=ai (even) 
U=8 (odd) 
U=oo (odd) 




FIG. 4. Local density {hi{t)) vs. time, showing local relaxation. 
Shown is the time evolution of an even and an odd site for [/ = 8 
and U — oo. The agreement is almost perfect. 
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FIG. 5. Time evolution of local densities on even sites for various U 
compared to the limiting behavior U — 0. 




1 2 3 4 ^5 



FIG. 6. Time evolution of local densities on even sites for various U 
compared to the limiting behavior U — oo. 



Under the assumption that for the time scales achievable by 
DMRG an asymptotic power-law can already be read off, one 
finds the exponents given in Fig. |7] While for very small U 
and all [/ > 4 one finds a slope similar to t^^^^, as for the lim- 
iting cases, there is an intermediate regime where relaxation 
is much faster. It must be stated, however, that in this regime 
the precise slope is hard to extract, even an exponentially fast 
decay cannot be excluded completely, but we are not aware of 
a physical reason why there should be a qualitative change of 
decay behavior from power-law to exponential. Qualitatively, 
stronger interactions should lead to a more distinct dispersion, 
leading in turn to a quicker decay. 

Concerning experimental implementations, a number of 
further issues require a discussion: The experimentally acces- 
sible systems will be of finite size, leading to recurrence ef- 
fects. Moreover, there will be a parabolic confining potential, 
which will modify results. Furthermore, one will encounter 
initial states at non-zero temperature, as well as possible ad- 
ditional fluctuations due to inhomogeneities in applied fields, 
different from the apparent local relaxation. 

So far, we have confronted analytical results in the i — > oo 



FIG. 7. Estimated negative exponents of asymptotic power-law de- 
cay law for local densities. The asymptotic decay exponent for the 
limiting cases is shown as a dashed line. 

limit with DMRG results for L = 32. In order to check 
whether on the time scales reachable by DMRG recutTence 
effects can be seen for this system size, we have rerun selected 
calculations for L = 50, observing no relative change in re- 
sults above 1 percent. In particular, the shapes of oscillatory 
behavior remained completely unchanged. In the U — limit, 
finite system DMRG results and infinite as well as finite sys- 
tem analytic results agree completely on the finite times reach- 
able by DMRG, even despite the different boundary condi- 
tions (open boundary conditions (OBC) for DMRG, periodic 
boundary conditions (PBC) for analytical statements). Recur- 
rences, where the difference between finite and infinite system 
becomes obvious, happen much later, see F ig.[8| 

The effect of the trap, as shown in Fig. |10) is to generate 
effective reflections from the edges of the system, leading to 
much earlier recutTences. However, in experimental setups, 
this effect would for realistic traps set in late enough to allow 
for sufficiently long observation. 

V. TIME EVOLUTION OF LOCAL CORRELATORS 

Let us now consider the nearest-neighbor correlator 
{ii+i(t)bi{t}) ■ This quantity is specifically interesting, and 
can in principle be measured by means of exploiting the tun- 
ing of the double well potential of the superlattice, and appro- 
priate timing (see above, and compare also Refs. |57-59|). It 
is of interest as it goes beyond local densities: The build-up of 
correlations in time starting from the uncorrected initial state 
becomes visible. 

In the limiting free cases, identical results are found: 

k=i (29) 

^-^^^Ji(4Jt) (L^oo). 

The real part of the correlator is strictly zero for all times, 
whereas the imaginary part relaxes to with an asymptotics 
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FIG. 8. Densities in the free limit: Analytics for L = 50 (PBC), 
L = oo (PBC) and DMRG for L = 50 (OBC; in fact calculated in 
the hardcore limit, because DMRG fails for U — 0, but densities are 
identical in both cases). On the time scales reached by DMRG, all 
three agree. This incidentally shows that boundary conditions are of 
little importance for the sizes and time scales of DMRG. Recurrence 
effects become visible later. 



FIG. 10. Density (ni(t)) (left) for [/ = 0, L = 100, includ- 
ing a trapping potential corresponding to a local chemical potential 
fii = 0.01 (j - L/2f. The plot on the right shows {Ne{t)) (black) 
and {No{t)) (green). Note the similarity to a Bessel function for 
short times and the recurrences after t « 10. The recurrences are 
even more pronounced in the quasi-momentum distribution, Fig.|18[ 
which lacks, however, signatures of local relaxation dynamics. 




FIG. 9. Density {n^{t)) (left) for [/ = 0, L = 100, and periodic 
boundary conditions. The plot on the right shows {Ne{t)) (black) 
and {No{t)) (green), compare Fig.[lo] where the same is shown in- 
cluding a trapping potential. 
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FIG. 11. Real part of correlations to neighbors as a 

function of time, for different values of U . Note that in the effectively 
free cases (7 = and U = oo, the absolute value of these correlators 
converges to zero. 



of t^^l"^ after a quick growth to a maximal value of about 
0.28 at time t « 1/2. This quick growth reflects the buildup 
of correlations due to particle motion with speed linear in J. 

For finite I] , the scenario has marked similarities and dif- 
ferences. Considering Figs.[TT|and[T2] one sees that on short 
time scales the buildup of the imaginary part of correlations 
is identical for arbitrary U . This simply reflects the fact that 
due to the distance 2 between particles at i = 0, no collisions 
have yet happened on these time scales. Only when the in- 
teraction strength becomes visible, the relaxation to follows 
different paths. As for local densities there is a clear trend that 
relaxation is fastest around [/ w 3, reflecting the particularly 
efficient scattering there. In the real part, convergence to a fi- 
nite value occurs for all finite [/, but not such a clear picture 
of relaxation speeds occurs. 



If one plots the asymptotic values of the real part of the 



nearest neighbor coiTelators (Fig. 13 i, one sees that there is 
a maximum around [/ « 3, reflecting the efficient scattering 
there. Interestingly, the large-?/ behavior can be very well 
approximated by a C/^^ curve for all ?7 = 4 and larger. 

Indeed, this dependence is exactly what one would expect 
in the thermal or Gibbs state of the Bose-Hubbard Hamilto- 
nian R: We hence look at local correlations in the Gibbs state 



-0H 



/Z 



with Z = tr[c"^^] . We will take 



i=l 



1). 



(30) 



(31) 
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FIG. 12. Imaginary part of correlations to neighbors {bl_^_-l(t)hi{t)) 
as a function of time, for different values of U. 



as our unperturbed Hamiltonian including interactions and 
look at leading orders in the hopping term 



(32) 



Then we find, using standard thermal perturbation theory, up 
to first order in J 



{bTb,)^kjm)- Idx 



Z 



(33) 



This means that we have ((5(i j) = 1 if i and j are nearest 
neighbors, zero otherwise) 



(34) 

g/3(7(n-m-l) _ ^ 

z^{n — m — 1) 



where z = q~^e„ -g^ _ Un{n — l)/2 are the local 
energies of the unperturbed Hamiltonian. So, indeed, within 
the validity of perturbation theory, we do find the anticipated 
linear dependence on 1 /[/, as seen also in DMRG simulations 
in the time-dependent scenario. By definition, the correlators 
merely probe local quantities. This corroborates the intuition 
that locally, the system is indistinguishable from the situation 
as if globally the system was in a state maximizing the en- 
tropy, under the constraints of motion |3j|5]0|. 

To complete the picture, let us finally remind ourselves of 
how the maximum entropy states of the global system locally 
look like. The constants of motion of interest here are the oc- 
cupation numbers in bosonic or fermionic momentum space, 
respectively. The global maximum entropy state — ^under the 
constraints of motion — for [/ = is then found to be 



Q = 



- In(3)f)j6i 



(35) 



This result is completely consistent with our result for the re- 



laxation of the characteristic function (see Appendix XIII B i 



FIG. 13. Equilibrated value of the real part of correlations to neigh- 
bors {bl^i{t)bi{t)) for large t for different values of U. As a guide 
to the eye, large-f/ behavior is fitted by a solid line proportional to 
1/U. 



One can now calculate expectation values in equilibrium for 
the local observables considered above (by construction in 
agreement with our earlier analytical findings) as 



(bib 



ij, (nifij) ~ {ni){nj) = -S^^j. (36) 



In the case of [/ — > cx) one finds 



which means that 



(37) 



(38) 



We have also numerically investigated the relaxation behav- 
ior of longer-ranged correlators (pj{t)bi{t)) . As an exam- 
ple, we depict in Fig. 14 the relaxation of the real part of 



(6|_i_2(i)6i(i)), the imaginary part being zero. Relaxation here 
is of monotonically increasing effectiveness with U . The 
relaxation dynamics for C/ = and C/ = oo is different, 
revealing the fundamentally different character of the non- 
interacting limiting cases. 



VI. LOCAL RELAXATION 

Let us summarize the properties of local relaxation that can 
be extracted in this experimental setup. In all quantities, there 
are three time regimes: 

1 . The first time regime is associated with the building up 
of correlations: This is governed directly by the cou- 
pling strength to the nearest neighbor, and on the short- 
est times is identical for all U, before collision pro- 
cesses become important. The time scale is Jt < 1, 
as this is when the first collisions happen and establish 
correlations. 
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FIG. 14. Real part of correlations to neighbors {bl^2(^)^iO')) ^ 
function of time, for different values of U. 



The second time regime is associated with local re- 
laxation. There is fast oscillating dynamics between 
neighboring sites, happening at a time scale dictated by 
the hopping J, which also governs the speed of sound. 
Local relaxation as such is slow, in the exact analysis 
a slow polynomial decay: The Bessel function fulfills 
I'^ols^)! < a:^^^^'' as a strict bound, with an asymptotic 
envelope 



-l/2^ 



(39) 



one finds relaxation {not due to decoherence but due the 
dilution of information over the lattice) to the true local 
equilibrium. For finite U, numerics is consistent with 
polynomial decay, which for intermediate U seems to 
be much/flifer, but still polynomial. 

The intuition is that this local relaxation is due to in- 
fluences of excitations travelling with the speed of in- 
formation propagation from farther and farther lattice 
sites, broadened by dispersion. This means that this 
time scale, even in the interacting case ofU > 0, should 
be defined by the speed of information propagation in 
the system. Note that to date, there is no rigorous bound 
known to the speed of information propagation in the 
fully interacting Bose-Hubbard model (this being a con- 
sequence of the lattice sites having an unbounded num- 
ber of particles, see also comments in Refs. |43-45|). 
On physical grounds, yet, it should be expected to be 
similar to J, i.e. "ballistic" transport at the high energies 
provided by the quench, slightly modified by U (possi- 
bly the resulting bounds can only be formulated for spe- 
cific initial states having small local particle number). 
This is also what the numerics shows. Note that the 
low-energy speed of sound need not be relevant here. 

The third (very large) time is the recurrence time, which 
seems to be shortened substantially in the presence of a 
realistic trap. As in the presence of a trap, the excita- 
tions no longer travel with a constant speed (the speed 
of sound), but are slowly reflected, one should expect 



a quicker relaxation, and this is also the behavior the 
analytics exhibits. Generally, for moderate to large sys- 
tem sizes it is already beyond the reach of our simula- 
tions. It would be interesting to see, possibly in exact 
diagonalization, whether the recurrences are weakened 
compared to the free solutions due to interactions in the 
system, see also Ref. fT3l. 

While these findings are essentially independent of the 
interaction strength U, we also find three local relaxation 
regimes for different U : 

• For very small values of U (up io U ~ 1) relaxation 
dynamics is quite close to the non-interacting bosonic 
limit. 

• For larger values of U the system seems to show a 
behavior very similar to the hard core bosonic or free 
fermionic limit with similar relaxation exponents of or- 
der 1/2. Local correlations are consistent with locally 
equilibrated subsystems perturbatively coupled in 1/U . 

• The case of intermediate U ^ i appears to be a special 
case for various observables, marking the "boundary" 
between the "free bosonic case" U — Q and the "hard 
core boson case" U — oo, so the fermionic one. Colli- 
sions lead to the most efficient relaxation in this case. 

To summarize, the dynamics of local quantities shown is con- 
sistent with the limiting U — Q and U — oo cases, but shows 
a richer phenomenology in particular for intermediate U . 



VII. QUASI-MOMENTUM DISTRIBUTION 

In this subsection, we briefly compare what is seen in the 
above local relaxation with the situation in momentum space. 
The quasi- momentum distribution. 



S{q,t) 



1 



E 



,i«(i-i)/!,t 



(40) 



obtainable from time-of-flight experiments, is no longer a lo- 
cal quantity, but a global one probing the state of the entire lat- 
tice, as well as the boundary conditions and possibly a confin- 
ing additional potential. In fact, in the free model for [/ = 0, 
one finds that the quasi-momentum distribution will be very 
little time-dependent, and if so, in a chaotic fashion showing 
little structure, see Fig. [17] More specifically, we find (see 
Appendix) S{q, t) = 1/2 for q = 2ttI/L, / = 1, . . . , L, and 



S{q,t) 



1 

L2 



L 



4 Jit cos(27rp/L) S,VP?{Lq/2) 



sin{2TTp/ L — q) 



(41) 



for all other q G [0, 27r]. In case of U = 0, the quasi- 
momentum distribution hence indeed does not show any sig- 
natures of the local dynamics, see Fig. 17 This is, again, not 
inconsistent at all with the notion of local relaxation: This 
quantity is a global one (and local only insofar as correlators 
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FIG. 15. Quasi-momentum distribution as a function of time, for 
U = 1.5. Tlie labels q on the x-axis correspond to momenta 
7r/(L + l)q, best suited to DMRG open boundary conditions, as ex- 
plained in the text. 
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FIG. 16. Quasi-momentum distribution as a function of time, for 
U = 3, marking the cross-over regime. Labels as in Fig.|15| 



of far away lattice sites are suppressed via a quickly rotating 
exponential function). For any subblock, the dynamics drives 
the system towards the values of equilibrium. In momentum 
space, this is not necessarily the case; at least not on these 
time scales. In the free situation of U = 0, this leads to an 
absence of relaxation of this quantity. 

It would still be interesting to measure this quantity for 
close to free cases for short times, to demonstrate the very 
point that the information about the initial condition is fully 
contained in the system, and merely diluted. Hence, obvi- 
ously, one cannot expect true global relaxation to happen, un- 
less a further mechanism of decoherence due to additional ex- 
ternal degrees of freedom is present. 

In Fig. 



18 we show the quasi-momentum distribution for 
the situation including a trapping potential. The effect of the 
additional harmonic potential is very significant. Yet, again, 
note the absence of signatures of the time scales of local dy- 
namics as compared to Figs. [9 10 



How does this compare to the case of having a finite on-site 
interaction C/? Again, the behavior reminds of the situation 




FIG. 17. Quasi-momentum distribution S{q, t) for U = 0, L = 
100, and open (left), periodic (right) boundary conditions. Note the 
dramatic difference between open/periodic boundary conditions and 
the absence of the time-scales clearly visible for local observables 
as in Figs. |9] 10 The fast oscillatory behavior and the relaxation 
time. See Fig. 18 for S{q, t) in the presence of a trap. The non-local 
quantity of the quasi-momentum distribution detects the boundary 
conditions early. 



in the free case of J7 = 0. The quantity is too non-local to 
grasp the effect of local relaxation. Still, for longer times, 
one may also expect relaxation in some translationally invari- 
ant properties in momentum space, and the quasi-momentum 
distribution may well relax. Such a situation is observed in 
a flow equations approach in case of an interacting model in 
Ref. ll30ll . For the DMRG calculations we use a slightly dif- 
ferent definition of the quasi-momentum distribution 



S{q,t) 



nqi 
L + l 



L + l 



(42) 

This ensures that S{q, t) is a constant of motion for q — 
1, . . . ,L in the free system ([/ = 0) in case of OBC. Fig. 
15 and Fig. 16 show this quasi-momentum distribution for 



q = 1, .... L for different values of U. Indeed, signatures of 



such a behavior are also observed in Fig. 15 where the mo- 
mentum distribution is depicted for different values of time t. 
For larger times, the quasi-momentum distribution appears to 
relax in the interacting case of U — 1.5. This effect is even 
more distinct in the case of stronger interactions, as depicted 



in Fig. 16 



To reiterate, the very absence of relaxation in close to free 
settings on short times gives rise to an interesting situation: 
One could experimentally see signatures of local relaxation. 
One could, however, also see that the information becomes 
more and more dilute in position space, but is still preserved 
in the system as such, as seen in the quasi-momentum distri- 
bution. 



Vni. GLOBAL DENSITY-DENSITY CORRELATOR 

One potentially experimentally accessible global correlator 
is the connected density-density correlator 



(43) 



which we show for different U in Fig. [19] Numerics indicates 
that this quantity looks effectively relaxed after some time as 
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IX. ENTANGLEMENT AND ENTROPY GROWTH 




FIG. 18. Quasi-momentum distribution S{q, t) for the same parame- 
ters as in Fig. [TO] Note tiie absence of the time-scales clearly visible 
for local observables as in Figs.|9][T0] The fast oscillatory behavior 
and the relaxation. The recurrence is however clearly visible. 
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FIG. 19. Global density-density correlations for a system of length 
L — 32 for various interaction strengths U. 



long-range correlators contributing will be small. For the free 
case of U = 0, one analytically finds for large L (see Ap- 
pendix) 



{N,{t)No{t))-{N,{t)){No{t)) 

L (44) 

^ - j^ (3 + M8Jt) - 4[Jo{4Jt)]^) , 

which relaxes to — 3L/16 for large times. For hardcore bosons 
the global density-density correlator for large L is given by 



{N,{t)Noit)) ~ {Ne{t)){Noit)) 
which relaxes to — L/16 for large times. 



(45) 



The creation of excitations at all points of the lattice will 
in general create a significant degree of entanglement in the 
time evolving state. Since at each site, an excitation starts 
to travel through the lattice, a linear increase of the entan- 
glement entropy of subblocks in time is to be expected ||2TI - 
[231 . More precisely, it has been shown in Refs. Il22l|23l that 
for any time evolution under a local Hamiltonian with finite- 
dimensional constituents, starting from a product state, the en- 
tanglement entropy of a subblock / ~ {1, . . . , s} grows in 
one-dimensional systems in s at most as 



(46) 



for suitable constants cq, ci. This means that for any constant 
time, the entanglement entropy satisfies what is called an area 
law. For larger and larger block sizes, the entanglement en- 
tropy will eventually saturate in s (and, e.g., not logarithmi- 
cally diverge). 

In turn, this upper bound is saturated, in the following 
sense: In Ref. 1 3 1 and explicitly in the Appendix 1-norm con- 
vergence is shown to products of Gaussian states 



||e/W-^max||l = ||wW-eil|l^O, 



(47) 



when considering initial conditions |-!/)(0)) ((0)| and time evo- 
lution under Hu=(), where pmax is a Gaussian product state. 
This observation immediately implies a bound to the entan- 
glement entropy that is linear in time, if the block size s can 
only be chosen appropriately. This has also been made ex- 
plicit in Ref. |33|, in that for the fermionic instance of this 
problem (or, equivalently, a spin model), there exist constants 
C2 , C3 , C4 , Lq , So , > such that 



Eigiit)) > C2t + C3, 



(48) 



for L > Lq and s > sq and to < t < C4S, again for 
I — {!,..., s}. In other words, for a larger and larger block 
size, one encounters a linear increase of the entanglement en- 
tropy of that block. In both Ref. |3| and |33| the local states 
converge to maximal entropy states under the constraints of 
motion, in the latter case starting from a fermionic Gaussian 
state. 

This means that eventually, one will have to use matrix- 
product states in the DMRG approach that make use of expo- 
nentially many parameters in time. This linear increase, prov- 
ably correct in the above cases, is consistent with a wealth of 
numerical findings in quenched systems, as well as with argu- 
ments using conformal field theory 1,15. 24 , 25 , 65 1 . 

Knowing that we have to asymptotically deal with a lin- 
ear increase, we have plotted the entanglement entropy as ob- 
tained from a t-DMRG approach. We depict here the entan- 
glement entropy as a function of time in a slightly different ge- 
ometrical setting, which still has implications on the approx- 
imatability of a state with a matrix-product state in DMRG. 
This is the setting of the symmetrically bisected half-chain, 
for s = L/2 and L — 32. The above linear lower bound 
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FIG. 20. Von-Neumann entropy of the bisected half-chain between 
sites {1, . . . , 16} and {17, . . . , 32} vs. time for various interaction 
strengths U . 



is also rigorously true for this geometry (at least for finite- 
dimensional constituents), whereas the linear upper bound is 
certainly expected to be valid. 

We numerically find an initial sublinear regime in time, fol- 
lowed by a linear regime, see Fig. 



20 The linear regime is 



plausible when considering the linear propagation of the exci- 
tations in the lattice. This sequence of a sublinear regime fol- 
lowed by a linear one is plausible when considering the obser- 
vation that eventually, excitations travel through the lattice at 
a finite speed. This is also consistent with the above linear up- 
per and lower bounds. This increase also eventually limits the 
time to which the time evolution can be traced using t-DMRG. 
It is interesting to observe that the entanglement growth de- 
pends only very weakly on interaction strength U, whereas 
the (low-energy) speed of sound depends quite strongly on U 
ll67l . This reflects the fact that after the quench, we are dealing 
with a very high energy state of the new Hamiltonian, where 
excitations are expected to be essentially ballistic with a prop- 
agation velocity cx J, such that the conventional low-energy 
speed of sound is not relevant. 



X. RELATIONSHIP TO KINEMATICAL APPROACHES 

We only very briefly touch this issue here. In this work, 
similarly to Refs. Il3]|5jpl, we have considered local relax- 
ation to an apparent equilibrium state. The overall system un- 
dergoes time evolution under a local Hamiltonian H, whereas 
subsystems / appear relaxed. This observation is in an in- 
teresting relationship to kinematical approaches 1371 147H491 . 
Indeed, if one draws a random pure state of a large system, 
this will be — in the limit of a large environment — almost al- 
ways maximally mixed. Random here means, drawn from the 
unitaiily invariant measure, so the Haar measure t47l450l (or a 
measure on the energy shell for Gaussian states |5T1). So one 
could argue that "most states are locally relaxed". This might 
render a relaxation to high entropy states plausible. Yet, ironi- 
cally, the image of the positive time axis under the time evolu- 



tion corresponds, of course, to a one-dimensional manifold in 
state space: So in order to show that relaxation follows from 
such a kinematical argument, one has to demonstrate that for a 
given local Hamiltonian, time evolution ensures that the state 
of the global system stays within the typical subspace. This is 
an interesting programme and the proof constitutes an inter- 
esting challenge in case of interacting systems, which has not 
yet been completed. In case of free systems as in Ref. |i3J, the 
findings may indeed be interpreted in this way. 



XI. SUMMARY 



In this work, we have introduced a setting in which local re- 
laxation in quantum many-body systems can be probed, with- 
out the need of actually addressing single sites: This is the 
setting of a Bose-Hubbard model, in which preparation and 
read-out can be done with period-2 translational invariance, 
as can be achieved by exploiting optical superlattices. We 
have approached this idea by deriving analytical expressions 
for the relevant quantities in case of the free limits ofU — Q 
and U = CO in the Bose-Hubbard model, as well as using a 
systematic t-DMRG approach in the time-dependent interact- 
ing case. For the [/ = case, we presented in detail a true 
convergence of subsystems to the maximum entropy state in 
1-norm. In several ways, the interacting setting reminded of 
the non-interacting case, certainly concerning the mechanism 
of relaxation. The time-scales are different, however, showing 
a faster relaxation compared to the inverse square-root depen- 
dence. Also, the dependence on the interaction strength re- 
flects the same dependence in the corresponding Gibbs state 
of the Hamiltonian the system is quenched to. 

In this way, signatures of local relaxation can be measured 
with present technology: Local densities are found to relax 
to those of a quasi-thermal state on well-defined time scales, 
often with a similar to inverse square-root time dependence. 
More sophisticated measurements would reveal correlators 
and density-density coiTelations, again showing local relax- 
ation in a characteristic fashion. Hence, by exploiting period- 
2 translational invariance, one has a tool at hand that effec- 
tively can be viewed as if one looked at a small local subsys- 
tem, showing all signatures of the local relaxation. In turn, 
the quasi-momentum distribution is a global quantity, one that 
shows in the free limit of U = no relaxation at all, and 
merely detects the boundary conditions quickly. 

Hence, this is a quite exciting situation that both the pres- 
ence of the memory of the initial condition could be experi- 
mentally probed — in the absence of relaxation in "too global 
quantities" — as well as local relaxation: Locally, the system 
"appears relaxed", for very long times, until recurrences be- 
come relevant. The technology offered by recent advances in 
the use of optical superlattices should hence open up the way 
to experimentally access such fundamental and old questions 
as the mechanism of local relaxation in quantum many-body 
systems. 
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XIII. APPENDIX: PROOFS 

We largely follow Ref |3| to show local relaxation. For 
the special situation at hand — one spatial dimension and 
[■(/"(O)) ~ |1, 0, . . . , 1, 0) — the proofs simplify significantly 
and we state them here for completeness. 



A. Preliminaries 

The bosonic operators evolve in time according to (this can 
be shown by solving Heisenberg's equation of motion or ap- 
plying the Baker-Hausdorff formula) 

L 

h{t) = e"^6,e-'*^ = ^ F,,,6,, V{L,t) = e'''" , (49) 
i=i 

where the entries of the Hamilton matrix H are given by 
Hi j = — (5dist(i For the translationally invariant setting 
one has 



1 ^ 



k=l 



where 

Xk = -2cos(27rfc/i) 
are the eigenvalues of V. 



(50) 



(51) 



of studying the quantum state, we investigate its characteristic 
function in phase space. Pointwise convergence of the char- 
acteristic function implies convergence in trace norm for the 
state. For subsets / C £ = {1, . . . , L}, the local state on / is 
given by a partial trace 



gj{L,t)^tic\i[\m){m\], 



(52) 



and the corresponding characteristic function to represent the 
state gj is defined as 



D, 



(53) 



Here the a.^ € C are the complex phase space coordinates and 



\m) 



-itH 



|1,0...,1,0) 



(54) 



is the time evolved state vector. We aim at showing that the 
state of any subblock / — {1, . . . , s} of s consecutive sites 
locally relaxes. For any such block, we find from Eq. ([49| that 



X=(l,0,...,l,0|nA(A)|l,0,...,l,0), (55) 
where we defined 

s 

A(i)=^a,y;,W. (56) 

We thus find the following explicit form of the characteristic 
function of t): 

L/2 L/2 

X = n<oi^2.(fe)io) n(ii^2.-i(/32,-i)ii) 



i=l 



1=1 



(57) 



where we made use of the explicit matrix elements of the dis- 
placement operator in the Fock basis. Now, /3^/3 = a.^a fol- 
lows from unitarity of V and we proceed by proving that the 
above product converges pointwise in a to exp(— q:^^q:/2). 



1. The causal cone 



B. Proof of local relaxation 

In this subsection, we rigorously show that for large times 
and large system size, the state of a subblock / becomes ex- 
actly a maximum entropy state. The key ingredients to this 
proof are Lieb-Robinson ideas and a central limit-type the- 
orem. The proof is not too technically involved, but also far 
from being trivial, and we present it for completeness. It is, af- 
ter all, quite astonishing that one does arrive at maximum en- 
tropy states without a time average. This proof largely follows 
Ref 1 3 1, adapted to our situation of an alternating sequence of 
bosons and no bosons per site as initial condition. 

We can now no longer merely think in terms of moments, 
as we want to prove convergence of the state itself Instead 



The key ingredient for most of what follows is the fact that 
the entries Vi^j — Vi^j become arbitrarily small for suffi- 
ciently large L and t. This can be seen as follows: The 
Vi,j = Vi-j may be thought of as Riemann-sum approxima- 
tion to the integral 




^^g2iJtCOs(0)gi/0 ^ i' J;(2Ji), (58) 



where Ji is a Bessel function of the first kind, for which one 
has Ji{x) < x~^^"^ for all a; > and all The error involved 
in this approximation is 

\V,-i^M2Jt)\<^-^^^, (59) 
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i.e., we have the bound 



(60) 



which converges to zero if we fix let L oo and then 
t — > oo. However, we need a bound on the entries of V for all 
I. To this end we complement the above bound with a Lieb- 
Robinson type bound on the influence of sites with large /: 
As Hi^j — for dist(z, j) > 1, we have {H^)i^j = for 
n < dist(i, j) =: d, i.e.. 



(61) 



Now, for any matrix M one has \Mi,j\ < \\M\\, where || • | 
indicates the operator norm. Furthermore, nl > (n/3)". We 
thus find 



\V^.J\ < 



< 



(6Jt)" _ {QJt/dY 
^ 1 - HJt/d 



^ d" 



, (62) 



independent of L. Hence, matrix entries Vij with dist(i, j) > 
6Jt are exponentially suppressed in dist(i,j), defining a 
"causal cone" as the influence of sites with dist(i,j) > 6Jt 
is exponentially small in dist(i, j). For given e > we have 
\Vij\ < e if dist(i, j) > B{t), where B{t) is given by the so- 
lution to (6 Ji/B(t))^(*) = e{l - 6Jt/B{t)). A crude bound 
on B{t) may be obtained from noting that 



{QJt/dY ^ QJt/d 



l-QJt/d - l-QJt/d' 



(63) 



Combining this with the bound for i,j inside the cone from 
above, we find for given e > and all i, j that 



where, using Eq. ( 50 1, 



r,s— 1 



L/2 
i=l 

(67) 



where we find for the last fine 

L/2 



L 
i.e., 

L/2 



J2 ^2.nis-r)/iL/2) ^ g^^^ ^ ^ 5._,,L/2, (68) 



-11! 

2L 



E' 

r=l 



,4Jit cos(27rr/L)g^ir(fc-;) 



(69) 



Thus, 



L/2 

|El/^: 



■2i| 



< J2 \ak\\ai\VkA2t)\, (70) 



\Vi A < e for ah i : — < Jt < 



47r2 8e + 6 ■ 



(64) 



The entries of V are thus arbitrarily small for sufficiently large 
L and t. In particular limt_>.oo limi_>.oo \ Vij \ = 0. 



which converges to zero for large L and t. We thus arrive at 
the desired statement: 



lim lim x(Q:;L,i) = c 

t-^oo L— foo 



(71) 



2. Central limit-type argument 



pointwise in a.. 



We have from the previous section that \Vij \ becomes ar- 
bitrarily small for sufficiently L and t. Thus, for given a. the 
< Sj=i I'^ill^'.il become arbitrarily small, in particu- 
lar we may assume \f3i\ < 1/2. Now, for |a;| < 1/2 one has 
I log(l + x) — x\ < x^, i.e.. 



L/2 

^|l0g(l-|/32._l| 



1/3: 



2i-l\ 



< 



i=l 



< 



L/2 

El/3: 

i=l 
L/2 

sup |/3,f ^1/32, 



2i-l I 



Now, 



L/2 s L/2 

E 1/^2,-11'= E alalJ2yk,2^-lVC2^-l, 
i=l k,l=l i=l 



(65) 



(66) 



C. Correlations after a sudden quench (U = 0) 

In this section, we give explicit forms for the two-point cor- 
relations. 



= {h\{t%{t)), 
and density-density correlations, 

gi,j = {n,{t)nj{t)), 



(72) 



(73) 



for the translationally invariant non-interacting case, U ~ 0. 
We already know the asymptotic behavior from the previous 
section (for long times, neighboring sites become uncorrected 
as the state converges to a direct product of Gaussians) but the 
goal here is to derive explicit expressions for all times. 
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We find from Eqs. ( |49|50| l that 



1 ^ 



^it{\,- \p) 2m{pj -qi) / L 



p,q=l 



L/2 



(74) 



27ri(2/c-l)(g-p)/L 



x-> e 

fe=i 



where the term in the last line evaluates to 



-27ri(9-p)/L 

L 



L/2 

fe=i 



27rifc(g-p)/(L/2) 



9 ~ "9-P,i/2 ~ 'Jp-q,L/2j 



(75) 



We thus have 



(-1)' 
2L 



p=i 



,4Jit cos(27rp/L) 27rip(i-i)/L 



(76) 



which yields a particularly simple form in the thermodynamic 
limit i — > 00 



1 



i-iy 



(77) 



where J„ is a Bessel function of the first kind. Thus, the total 
number of particles at even sites 



(78) 



is a truly local quantity in this translationally invariant setting. 

The above allows us to derive the quasi-momentum distri- 
bution. 



1 1 

2 " 2L2 



E< 



E< 

p=i 



,4 Jit cos(27rp/L) 



(79) 



jii(ij+7r-27rp/L) ij(2Trp/L-(j) 



where the last term vanishes for Lq/ (2tt) E {1, . . . , L}, yield- 
ing S{q, t) = 1/2. For all other q G [0, 2tt] we find 



p=l 



^4 Jit cos(27rp/L) 

sin2(Lq/2) 
sin(27rp/_L — (7) 



(80) 



To derive the density-density correlations, we first observe 
that for 

hkj,r,s = ("-1,^2, . . . ,nL\b\biblbs\ni,n2, ■ ■ .,nL), (81) 
we have (writing Sk^i = 1 - 4j) 

hk,l.r,s = {Sk,l + Sk^l) {Sr,s + Sr^s) hk,l,r,s 

= SkdSr^snkUr + 5k^l5r^shk,Lr,s (82) 

= h,iSr,snknr + 5k^iSk^sSi,rnki'^ + ni). 



Thus, 



L/2 L/2 
9i,] = E ^i*2k-l^i,2k~l E ^j*2r-l^",2r-l 



k=l 



L/2 



2 ^ 45^/^^i*2fe-l^i,2i-lV}*2;_iV,-,2fe-l, 



fe,/=l 



and therefore. 



(83) 



(84) 



(85) 



9i,j —fi,ifjj + fi,j{Si,j + fj.i) 
L/2 

-2^iK2fc-inis-,2fe-ir 

fc=l 

This expression yields 

L/2 

— (ff2i,2j-l ^ f2i,2if2j-l,2j-l) 

L/2 

= /2i,2j-l/2j"-l,2z 
L/2 L/2 

— E l^2i,2fc-inV2j-l,2fe-lP, 

k—1 — l 

where, using the unitarity of V, 

L/2 

E l^2.,2fc-inV^2,-l ,2fc-l| 

/ L/2 \ L/2 

= I l-El^2»-l,2fc-in El^2,-l,2fe-l| 

= (1 — /2fc-l,2/£-l)/2/£-l,2/c-l, 



and, using the explicit form of fij and identifying Kjonecker 
deltas, 

L/2 

f2i,2j-lf2j-l,2i 
^ 4Jit(cos(27rp/L)+cos(27rg/L))-27ri(p-g)/L 



(86) 



16 



P,9=l 



1 

8 Jit cos(27rp/L) 



16 



p=i 



L 
16' 



(87) 



For large L we then find 



L (88) 

-^^ (3 + M8Jt) - 4[Jo(4Jt)]') . 
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D. Correlations after a sudden quench (U oo) 



Using the results for U — (XIIICi we can now easily 
derive analytical expressions for local densities 



next-neighbour two-point correlations 

and density-density correlations 

9i,j = {n.,{t)nj{t)) 



(89) 



(90) 



for [/ — > oo. For simplicity we assume even particle numbers 
N = L/2. With Eqs. (|22|, ^ it follows that Eqs. ([89 1- 
( [9T] i are also valid for after the mapping to fermions, with 
the bosonic creation- and annihilation operators replaced by 
fermionic ones. Because of the identical time evolution of the 



bosonic- and fermionic operators (equation ( 12 1 respectively 
( [25] l) we find essentially the same results for the local densities 
and the next-neighbour two-point correlations like in section 



XIII C as we did not make any use of commutation relations 



in the derivation of the bosonic respectively fermionic results: 



J^,^ 2 2L 



li+lA 



i-ir 

2L 



L 



,4Jit cos(27rp/L) 



,4Jit cos{2iTp/L) -2mp/L 



(92) 



(93) 



In contrast to this, the results for the density-density correla- 
tions differ For fermionic operators a similar derivation to Eq. 



( 82 1 yields 



. ,nL\clciclcs\ni,n2, ■ 



,nL) 



This finally leads to 

9i,j — fisfj.J + fi.j{Si,j — fj,i) 



(94) 



(95) 



with fij = {cl{t)cj{t)). Using the explicit form of the time- 
evolved operators in Eqs. ( [T2] i, (25 i and the absence of any 
signs resulting from commutation relations it i s easv to see 



(91) that fij is identical with fij for free bosons (see 



finally leads to Eq. (26i. Using the results of XIII C the global 



XIIICi. This 



density-density correlator is then given by 

L/2 

L 

8Jitcos(27rp/L) _ 



(96) 




For large L this yields 



L 
16 



(1- Jo(8Ji))- 



(97) 
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